Momentum-resolved radio-frequency spectroscopy of ultracold atomic Fermi gases in 

a spin-orbit coupled lattice 
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We investigate theoretically momentum-resolved radio-frequency (rf) spectroscopy of a non- 
interacting atomic Fermi gas in a spin-orbit coupled lattice. This lattice configuration has been 
recently created at MIT [Cheuk et ai, arXiv: 1205.3483 for 6 Li atoms, by coupling the two hy- 
perfine spin-states with a pair of Raman laser beams and additional rf coupling. Here, we show 
that momentum-resolved rf spectroscopy can measure single-particle energies and eigenstates and 
therefore resolve the band structure of the spin-orbit coupled lattice. In our calculations, we take 
into account the effects of temperatures and harmonic traps. Our predictions are to be confronted 
with future experiments on spin-orbit coupled Fermi gases of 40 K atoms in a lattice potential. 
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I. INTRODUCTION 

The past few years have witnessed an exponential 
growth of interest in studying ultracold atomic gases 
under a synthetic gauge field [H4l6|. The growth is 
strongly motived by a series of ground-breaking exper- 
iments at the National Institute of Standards and Tech- 
nology (NIST) [H-[i|. Most notably, synthetic spin-orbit 
coupling - the coupling between the spin and orbital de- 
grees of freedom of the atom - was created and detected 
in an atomic Bose-Einstein condensate (BEC) of 87 Rb 
atoms in early 2011 [2j. Such a spin-orbit coupling is re- 
sponsible for the recently discovered topological states of 
matter, such as topological insulators and spin quantum 
Hall materials [t7| . [l8| which are new types of functional 
materials that may lead to novel quantum devices. It is 
natural to anticipate that the investigation of spin-orbit 
coupled ultracold atomic gases will provide an entirely 
new platform to simulate and understand new genera- 
tion materials. 

To date, spin-orbit coupled atomic Fermi gases have 
been realized at ShanXi University [6( and at the Mas- 
sachusetts Institute of Technology (MIT) [7], by using 
fcrmionic 40 K atoms and 6 Li atoms, respectively. The 
technique used to induce spin-orbit coupling in Fermi 
gases is more or less the same as in BECs. A pair of 
counter-propagating laser beams along x-axis is used to 
connect two atomic hyperfine spin-states, labeled by |t) 
and via a two-photon Raman transition. The Ra- 
man beams impart momentum 2hk^e x to a fermion while 
changing its spin from \l) to |^)- In this way, the orbital 
motion is coupled to spin and an effective spin-orbit cou- 
pling is generated. In the MIT experiment @, an ad- 
ditional radio-frequency (rf) coupling is used to couple 
the two hyperfine spin-states. Combined with the pair 
of Raman beams, this creates a periodic lattice poten- 
tial, in addition to spin-orbit coupling. The rich band 
structure of such a novel spin-orbit coupled lattice has 
been characterized through spin-injection spectroscopy 



, which uses a rf laser beam to inject free atoms in a 
third spin state into an empty spin-orbit coupled system, 
and then obtains the momentum and spin of injected 
atoms using time of flight and spin-resolved detection. 
The spin-injection technique is particularly useful for 6 Li 
atoms. Due to the rapid heating from Raman process, 
the spin-orbit coupled Fermi gas of 6 Li atoms can hardly 
be created in equilibrium. The heating problem can be 
avoided by the spin-injection of atoms from a free Fermi 
gas in the third spin state, which do not experience the 
Raman process. We note that for 40 K atoms the heating 
issue due to Raman process is much milder. As a result, 
a spin-orbit coupled Fermi gas of 40 K atoms can be cre- 
ated in equilibrium at T ~ 0.6Tp where Tp is the 
Fermi temperature. 

In this paper, we investigate theoretically momentum- 
resolved rf spectroscopy of a non-interacting, trapped 
atomic Fermi gas of 40 K atoms in a spin-orbit coupled lat- 
tice, given the perspective that such a system can easily 
be realized at ShanXi University [6] . This can be viewed 
as the first step to understand momentum-resolved rf 
spectroscopy of a strongly interacting atomic Fermi gas 
in spin-orbit coupled lattice. The momentum-resolved rf 
spectroscopy, whose initial state is a spin-orbit coupled 
Fermi gas in equilibrium, yields equivalent information to 
spin-injection spectroscopy for a non-interacting system. 
However, the latter approach may hardly give useful in- 
formation for a strongly-interacting system, due to the 
lack of equilibrium in the final spin-orbit coupled state. 
In our calculations, we take into account the effect of har- 
monic traps by using local density approximation. The 
effect of temperatures is also addressed. 

Our paper is organized as follows. In the next section 
(Sec. II), we give the model Hamiltonian of spin-orbit lat- 
tice and explain how to calculate the single-particle ener- 
gies and eigenstates. In Sec. Ill, we derive the expression 
for momentum-resolved rf-spectroscopy and discuss the 
results for a homogeneous spin-orbit coupled system. In 
Sec. IV, we present the rf spectroscopy of a trapped sys- 
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tern within local density approximation. We discuss in 
detail the evolution of the rf spectroscopy as functions of 
the temperature, Raman coupling and rf coupling. Fi- 
nally, Sec. V is devoted to conclusions and some final 
remarks. 

II. MODEL HAMILTONIAN 



and the rf lattice potential 

V L (x) = Vl [cos (2k R x) a x + sin (2k R x) a y ] . (11) 

Here, for convenience we have introduced a spin-orbit 
coupling constant A = h 2 kji/M, an "effective" Zeeman 
field h = VLr/2, and an "effective" lattice depth Vl = 
Qrf/2. 



A non-interacting atomic Fermi gas in a spin-orbit cou- 
led lattice may be described by the model Hamiltonian 
, H = Ho + Hr + Hrf, where, 



n = Yl J (r) 



H R = ^ / ,/r 



A. Single-particle solution for Hso 

The model Hamiltonian Hso describes a spin-orbit 
coupling with equal Rashba and Dresselhaus strengths 
0)@-0|- The single-particle solution 0k (r) satisfies the 
Schrodinger equation, i?so0k(r) =£k0k(r). Using the 
Pauli matrices and the fact that the wave- vector or mo- 
if)]. (r) e l ™ RX ipi (r) + H.c. , (2) mentum k = (k x ,k±) = (k x ,k y ,k z ) is a good quantum 



h 2 k 2 



2M 



Vv (r) 



(1) 



Urf — 



RF 



dr 



ip\ (r) ip± (i 



H.c. 



(3) 



Here, ipl (r) is the creation field operator for atoms in 
the spin-state a = |t) and ||). The Hamiltonians Hr 
and %rf describe, respectively, the pair of counter- 
propagating Raman laser beams and the additional rf 
coupling that couple the two hyperfine spin-states. £Ir is 
the Raman coupling strength, kR = 2tt/X is determined 
by the wave length A of two lasers and 2hkR is the mo- 
mentum transfer during the two-photon Raman process, 
Qrf is the rf coupling strength. 

The Hamiltonian Hr creates the spin-orbit coupling, 
and with Hrf, a spin-orbit coupled system in a lattice 
potential can be formed. To see this, let us take the 
following gauge transformation, 



^ t (r) 



_ „+ik R x 



0) > 



^(r) = e- 4fe ^(r), 
with which the model Hamiltonians become, 



Ho 

Hr 
Hrf 



h 2 (k±k R e x f 
2M 



^ (r) 
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dr 
dr 



$(r)^. (r) + H, 



(4) 
(5) 



,(6) 
(7) 



$ (r) e -«**»^. (r) + H.c.l , (8) 



where in the first term of Hq we take "+" for spin- up 
atoms and "— " for spin-down atoms. By introducing a 
spinor field operator $(r) = (r) ,tp± (r)] T and using 
the Pauli matrices a xi a y , and a z , we can write compactly 
the model Hamiltonian in the form, 



H = J dr$t (r) [H so + V L (z)] * (r) , 

where we have defined the spin-orbit Hamiltonian 

k 2 ^ 



H 



so 



h 2 (4 



2M 



ha x + \k x <r z 



(9) 
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number, it is easy to see that we have two eigenvalues 
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-±Vh- 



X 2 k 2 



(12) 



2M 2M 
where "±" stands for two helicity branches. The cor- 
responding eigenstates are given by (we set the volume 

v = i), 



COS 0k \ ik x ., 

sin 0k 



4 -) ( r ) = 



- sin 0k \ e ik m x 



where 0k = arctan[(-\//v 



ik_L-r_L 



(13) 



COS 0k 

-' J- \ 2 k 2 - Xk x )/h] and rj_ 



(14) 



B. Single-particle solution for the spin-orbit 
coupled lattice 

In the presence of the additional rf Hamiltonian Hrf, 
the momentum along the x-axis, k x , is no longer a good 
quantum number. The lattice potential terms cos (2kRx) 

and sin (2kRx) will couple the eigenstates (r) and 

(r) if k' x - k% = 2nk R , where n = ±1,±2,--- is 
an integer. In this case, it is useful to define a quasi- 
momentum or lattice momentum q x for arbitrary k x as 
follows: k x = 2nkR + q x , where the integer n is chosen to 
make —kR<q x <kR. The quasi-momentum q x is then 
a good quantum number. We may expand the single- 
particle eigenstate of the total Hamiltonian in the form, 



-+-QO 



r) 



where k m = k^ + (2mkR 
same quasi-momentum q 
and 



Q_x)^x — kj_ -\- k mx e x 
and the energies of <ft£~' (r) 



(15) 
has the 

(+) 



k (r) are given by 



h 2 k\ 

2M 



R + kmx) 



2M 



±^h 2 + \ 2 k 2 mx . (16) 
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The coefficients a n + and a„_ can be determined by the 
Schrodinger equation, 

[Hso + V L (x)}$(q x ,k±;r) = E(q x ,k±)$ (q x , k ± ;r) , 

(17) 

where E(q x ,k±) = E(q x ) + h 2 k 2 L /{2M). By substituting 
the wave-function (|15p into the above Schrodinger equa- 
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tion and multiplying on both sides (fr^ ' (r) and finally 
taking the integration f dr, it is straightforward to show 
that, 



( e„+ \ / a n+ \ ^ / V++ V+ m \ ( a m+ \ 
{ e„_ J { a„_ J + J^^ \ V~+ V-- )\a m .) 

where the matrix elements of the rf Hamiltonian are given by, 

Vnm = v l (+<Wi,m cos 9 kn sin 6> km + 5„, m+ i sin 6» k „ cos 6> km ) , (19) 

Vnm = v l (+<Wi,m cos 6» k?i cos 6> km - (5 n , m+ i sin 6> kji sin 6> km ) , (20) 

V~+ = V L (S n +i,m sin 6» kn sin 6» km + J n , TO+ i cos 6» k „ cos 6» km ) , (21) 

Vnm = (-<5n+i,mSin0 k „ cosfl km - <5n,m+i cos6» k?i sin6» km ) . (22) 



2 1.^ 



2M 



+ £ (<?*) 



8n+ 
tin- 



(18) 



For a given quasi-momentum g^, therefore the wave- 
functions a n+ and a„_ and the corresponding energy 
E (q x ) can be obtained from the secular equation (fT5)) 
by exact diagonalization. In the numerical calculation, 
we have to make a truncation, i.e., let — ./V max < n,m < 
JV mia , where -/V max is a large integer. We will take and 
En = h 2 k 2 R / (2M) as the units for (quasi)-momentum 
and energy, respectively. 




Figure 1: (color online) Energy band structure E (q x ) in the 
absence (a) and presence (b) of spin-orbit lattice potential. 



C. Band structure of spin-orbit coupled lattice 

In Fig. 1, we present the band structure E (q x ) of the 
spin-orbit coupled system at Qrf — and £Irf — Er . In 
the absence of rf-coupling (Fig. la), the band structure is 
actually exactly identical to the single-particle dispersion 



Eq. (fT2"|) . However, we have folded the entire dispersion 
into the first Brillouin zone —kn < q x < kn. Thus, at 
the edge of Brillouin zone, q x = ±fc/j, the band energy is 
at least two- fold degenerate. Moreover, some bands are 
also degenerate at q x = due to the even-parity of the 
single-particle dispersion (fT2"j) . In this case, for the wave 
function Eq. (|15|) . there is only one non-zero (i.e., unity) 
coefficient in a„+ or a„_. With the rf-coupling (Fig. lb), 
the degeneracy at the zone edge or at q x = is lifted. 
We have a clear, well- resolved band structure. When 
the rf-coupling £Irf is large enough, a band gap also 
opens among the lowest three bands. Therefore, when 
the system is filled up to the top of the second band 
(i.e., the chemical potential fi < Er), a band insulator is 
formed. 



III. RADIO-FREQUENCY SPECTROSCOPY IN 
FREE SPACE 



Let us consider the rf-spectroscopy [19j, |20j, which is 
driven by a rf laser beam to transfer an atom in one of 
the two hyperfine states (say |!)) to an empty hyperfine 
state 1 3). The state |3) is normally higher in energy by 
an amount of hw^, due to the magnetic field splitting 
in bare atomic hyperfine levels. The Hamiltonian for the 
rf-transition may be written as, 



V rf = V I dr [4 (r) ^ (r) + (r) yj 3 (r)j , (23) 
= V [ dr \e- lk * x i:l (r) ^ (r) + H.c.l , (24) 



where ^3 (r) is the field operator which creates an atom 
in |3) at the position r and Vq is the strength of the rf 
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drive. In the second line of the above equation, we have 
taken the gauge transformation for ijj^ (r). As a result, 
there is an effective momentum transfer k R e x . 



The transfer strength of the rf-transition T (w) can be 
calculated by using the Fermi's golden rule: 



r ( w ) = J2 K $ /l v rf l^)l 2 / ( E i -ii)S[nu- nu H - (e s - ^)] . 



(25) 



Here, the summation is over all the possible initial states $i (with energy Ei) and final states (with energy Et) and 
/ (Ei — /j,) is the Fermi distribution function. The Dirac delta function ensures energy conservation during transition. 
Hereafter, without any confusion we shall ignore the energy splitting in the bare atomic hyperfine levels and set 
w&j, = 0. To calculate the overlap between the initial and final wave-functions |($/| V r f let us take the l-th band 

eigenstate 5>i = (q x , k^; r) = [>ipi' (r) , ijj^ (r)] T as the initial state, where 



^ (r) = ^2 \ a n+ sin ^k„ + O-n- cos 



e i(2nk R +q 3: )x £ i'k ± -r ± 



(26) 



n— — oo 



and Ei — E^ l \q x ) + h 2 k 2 _/(2M). It is easy to see that, in order to have a nonzero overlap, the final state must be a 
plane wave, i.e., e l (^-k a e x )-r ^ w j£ n wn ich the overlap of wave functions is given by, 



and the final state energy is 



<$/|Vr/|$i 



a^!l_ sin #k„ + ffli _ cos 9^ 



(27) 



h 2 k\ h 2 (k nx ~k R y 



2M 2M 

By taking into account all the possibilities for and the transfer strength can then be written in the form, 

oo + kj. 



(28) 



/, „ „ oo +oo 

{27T) { Hn=- 
-k R 

+ E^(q x )-fl 



-fcj 



a^l. cos 6> k „ 



x/ 



2Af 



w 2Af 



(29) 



where fc na: = 2nkn + q x and 6*k„ = arctan[(-\//i 2 + A 2 ^^ — Afc na: ) //i] . The integration over k± can be done analytically. 
We find that. 



k±dk± 



ft 2 *; 2 

2M 



+ £ (0 (<Z*)- 



MksT . 
4n 2 h 2 



In < 1 + cxp 



E^(q x 



k B T 



(30) 



Experimentally, the momentum of the transferred atom, k nx , could be resolved fj\, [2l|. Therefore, we may define a 
momentum- resolved transfer strength (k x = k nx — ks), 



T(k x ,uj) 



Mk B T 

47T 2 ft 2 



j^a„+ sm ^k„ + cos #k„ 



In < 1 + cxp 



E {l) {q x ) - /i 
k B T 



hLU + E^(q x ) ~ 



h 2 kl 



2M 



, (31) 



where the quasi-momentum and the index n are now 
determined from the momentum k nx = k^ + k R . The to- 
tal transfer strength is simply T(uj) — J_ dk x T(k x ,oj). 
To take into account the energy resolution of the spec- 



troscopy 7 ~ O.IEr (tJ, we may replace the Dirac delta 
function by, S (x) = (j/ir)/[x 2 + j 2 ]. To reveal more 
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clearly the band structure, it is also useful to calculate 

T(k nx ,u))=v[k x + k R ,uj + ^\ , (32) 
for which, the Dirac delta function takes the form S[huj + 

A. Momentum-resolved rf spectroscopy in free 
space 

We are now ready to calculate the momentum-resolved 
rf spectroscopy of a uniform Fermi gas in the spin-orbit 
coupled lattice. For a given k x , we obtain first the quasi- 
momentum q x and the index n, and then solve the energy 
bands E^ l \q x ) and eigenstates a^l and a^_. For the 
given chemical potential \i and temperature fc^T, which 
are in units of Er, we finally take the summation over the 
band index I and obtain the momentum-resolved transfer 
strength. Figs. 2a-2c and 2d-2f report the linear contour 
plot of rf spectroscopy in the absence {Qrf = 0) and 
presence (Qrf = Er) of lattice potential, respectively, 
with increasing the chemical potential fi at zero temper- 
ature. The transfer strength always becomes stronger 
with increasing fi, since there are more and more atoms 
in the system. 

The spectroscopy is already very interesting without 
a lattice, as shown in Figs. 2a-2c. At a low chemical 
potential \i — Er, only the lower helicity branch of the 
single-particle dispersion Eq. (|12[) is occupied. The spin- 
orbit coupling leads to a long tail at negative momentum 
and high frequency, in sharp contrast to a single Dirac 
delta function 5(oj) in the absence of spin-orbit coupling. 
As the chemical potential increases (Figs. 2b and 2c), the 
upper helicity branch gets occupied. In the spectroscopy, 
this creates a strong response at opposite momentum and 
frequency. 

In the presence of a lattice potential induced by the 
rf-coupling, the spectroscopy is greatly modified by the 
formation of energy bands. At the chemical potential 
fi = Er the lowest two bands, well-separated in energy, 
should already be occupied, as we can see from Fig. lb. 
As a result, we observe in Fig. 2d the two different re- 
sponses from the two bands. The energy gap between the 
first and second bands, which is at about Er, is clearly 
resolved in the spectroscopy. With increasing the chemi- 
cal potential, more and more energy bands come to con- 
tribute and the spectroscopy becomes more fragmented. 
There are some discontinuity at k x — ik R , indicating 
the existence of different energy gaps. 

Fig. 3 shows T(k nx ,ui) at /i = 5Er with or without the 
lattice potential. As T{k nx ,Cj) contains the Dirac delta 
function 6[hu) + E^(q x )], we anticipate that the energy 
band E®(q x ) can be seen clearly from the contour plot. 
Indeed, we find that the rf-response is peaked exactly at 
—By) (q x ) (shown by white curves) within the experimen- 
tal energy resolution. The strength of the response is de- 
termined by the coefficients aS_ and . Therefore, by 
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Figure 2: (color online) Momentum-resolved rf spectroscopy 
of a uniform Fermi gas without (left panel, Qrf = 0) or with 
(right panel, Q,rf = Er) lattice potential at zero tempera- 
ture. The Raman coupling strength is £Ir — 2Er. From (a) 
to (c), or from (d) to (f), we increase the chemical potential 
fi from Er to 5Er. Here the intensity of the contour plot 
increases linearly from (blue) to 0.1M/ft 2 (red). 
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Figure 3: (color online) Linear contour plot of T(k nx ,Cj) at 
Q R = 2E R , (j, = 5E R , and T = 0. The left plot and right plot 
correspond to the cases without and with lattice potential, 
respectively. The energy bands ~E^ l \q x ) are shown by thick 
white curves. 



measuring momentum-resolved rf-spectroscopy T(k x ,u>) 
and re-constructing T(k nx ,Q), we are able to obtain the 
complete information of the single-particle energy bands 
and eigenstates. 
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IV. RADIO-FREQUENCY SPECTROSCOPY IN 
HARMONIC TRAPS 

We turn to address the realistic issue of harmonic 
traps, Vr(r) = Mu) 2 r 2 /2, by using local density ap- 
proximation. Within local density approximation, the 
whole system can be regarded as a collection of many 
uniform blocks with a local chemical potential (i— Vy(r). 



The momentum-resolved rf spectroscopy is a sum of local 
spectroscopy over the whole trap, 

T T (k x , u) = J drT [k x , u;n- V T {r)] . (33) 

By substituting Eq. (|31|) for local spectroscopy and using 
the fact that 



Attt 2 dr\n < 1 + exp 



2nk B T 
Muj 2 



3/2 



5/2 -exp 



£ (0 (<?*)-/i" 



kuT 



we find that, 



(34) 



r T (k x ,u) 



E R (k B T) 



5/2 



7T (fiw ) d 



oo 



a«sin0 k co S 0, 



5/2 -exp 



2M 



(35) 



Here Li ra (cc) is the polylogarithm function. Experimen- 
tally, the trap frequency is about two order smaller 
than the recoil energy Er. Hereafter, we shall take 
Hloq = 0.01-Er, according to the experimental setup at 
ShanXi University [6j. In analogy to the uniform case, 
we may define 

TT(k nx ,u) = r T [k* + k R ,u + —z), (36) 



in order to better visualize the energy band. 
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Figure 4: (color online) Comparison between momentum- 
resolved rf-spectroscopy in free space and in harmonic traps 
at zero temperature. Here, we take the same chemical po- 
tential jj, — 3Er and use Qr — £Irf = 2Er. The intensity 
of each contour plot increases from (blue) to its maximum 
value (red). 

In Figs. 4a and 4b, we compare the momentum- 
resolved rf spectroscopy in free space and in harmonic 
traps, at Q,r = Qrf — 2Er and at zero tempera- 
ture. The chemical potential is taken the same value, 
i.e., [a = 3Er. In each plot, we assign the red color to 



the maximum value of rf transfer strength. With traps, 
the rf-response from the higher bands is blurred by the 
trap average. However, the qualitative features of rf- 
spectroscopy remains the same, as we may anticipate. 
This strongly indicates that in harmonic traps we could 
still be able to re-construct the energy band structure by 
using momentum-resolved rf-spectroscopy. 
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Figure 5: (color online) Total number of atoms as a function 
of chemical potential for a trapped Fermi gas. The cases with 
or without lattice potential are plotted by solid and dashed 
lines, respectively. For comparison, we show also the result 
without spin-orbit coupling by a thin solid line. Here, we take 
the trapping frequency fruJo = O.OIEr. 

To have a realistic estimate of the chemical potential in 
harmonic traps, it is useful to calculate the total number 
of atoms, which is given by N = J drn(r), where the 
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local density 



n(r) 



k±dkj 
o -k 



+k R 



1=0 



h 2 k\ 

2M 



By integrating over k± and the spatial coordinates, we obtain, 



N 



Er (k B T) 



5/2 



/ dfe^Li 5 



+ £«(<Z ;c ) -/i + V T (r) 



£ (i) (<?*)-M 



(37) 



5/2 - exp 



A: B T 



(38) 



In Fig. 5, we plot the total number of atoms in harmonic 
traps as a function of chemical potential with or with- 
out the spin-orbit coupled lattice. For comparison, we 
show also the result for an ideal Fermi gas without any 
spin-orbit coupling, N = /i 3 /[3(huo) 3 ]. In the ShanXi 
experiment [f|, the number of atoms is about 2 x 10 6 , 
corresponding to /i ~ 3Er. 



Momentum-resolved rf spectroscopy in 
harmonic traps 



We examine first how the rf-spectroscopy is affected 
by temperature. In Fig. 6, we report the evolution of 
rf-spectroscopy with increasing temperature ksT from 
0.2/i, 0.6[i to fi, where /i — 3Er, £Ir = 2Er, and 
Qrf — Er. As the temperature increases, more and 
more energy bands are visible, as these bands become 
thermally occupied. Importantly, there is no significant 
thermal broadening for higher energy bands. They are 
all well-resolved even close to the degeneracy tempera- 
ture ksT ~ (i. Note that, the typical temperature in the 
ShanXi experiment is about O.GfceTF ~ 0.6/1. 

We now explore the spin-orbit coupled system for a 
range of coupling strengths. In Fig. 7, we report the 
evolution of rf-spectroscopy as functions of the Raman 
coupling SIr and the rf coupling Qrf, at ksT = 0.6/x and 
fx = 3Er . This may be viewed as a realistic simulation of 
a future experiment for a non-interacting trapped Fermi 
gas of 40 K atoms in spin-orbit coupled lattice [6j . 

The reconstructed plot of TT{k nx , w) is shown in Fig. 
8, for a set of parameters £Ir — FIrf — 2Er, fi = 3Er, 
and T = 0.6/1. By comparing with the energy bands 
—E"\q x ), which is plotted by thick white curves, it 
is readily seen that the band structure can be clearly 
extracted from a realistic momentum-resolved measure- 
ment at finite temperatures (i.e., T ~ 0.6Tf) and finite 
energy resolutions (i.e., 7 ~ O.IEr). 




klk 



Figure 6: (color online) Temperature dependence of the 
momentum-resolved rf spectroscopy of a trapped atomic 
Fermi gas at Q.r — 2Er, Q,rf = Er and /1 = 3Er. The 
intensity of the contour plots increases from (blue) to 
2.5S|/(^ ) 3 (red). 



V. CONCLUSIONS 

In summary, we have predicted theoretically 
momentum-resolved rf spectroscopy of a non-interacting 
atomic Fermi gas in a spin-orbit coupled lattice. We 
have shown that such a rf-spectroscopy, just like the 
angle-resolved photoemission spectroscopy in condensed 





-3-2-10123 



k Ik 
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Figure 8: (color online) Linear contour plot of rT(fe n x,w) 
at Q, R = n RF = 2E R , n = 3E R , and T = 0.6/Lt. The 
energy bands —E^ l \q x ) are shown by thick white curves. 
The intensity of the contour plots increases from (blue) to 
l.<oE R /(hw f (red). 



matter physics [22j, provides an ideal technique to 
characterize the non-trivial band structure of spin-orbit 
coupled lattice. Our predictions can be readily examined 
at ShanXi University by using ultracold 40 K atoms [6|]. 

Using Feshbach resonances [23j], a strongly interact- 
ing Fermi gas of 40 K atoms in spin-orbit coupled lat- 
tice would be created experimentally very soon. It is 
of great interest to study how the single-particle band 
structure is modified by strong interatomic interactions 
and fermionic superfluidity. In this case, we anticipate 
that momentum-resolved rf spectroscopy would provide 
very useful information. 
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